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ABSTRACT 


This  is  a user's  manual  for  the  TIMES  computer  program, 
an  integrated  system  of  FORTRAN  subroutines  for  modelling  of 
a single  time  series  and  for  transfer  function  modelling 
with  a single  input  series.  TIMES  is  an  extension  and  adap- 
tation of  a widely  used  system  available  from  the  Ohio  State 
University  Computer  Center  Library.  This  document  is  a re- 
finement of  an  earlier  guide  prepared  for  classroom  use 
during  Spring  Quarter,  1976,  at  the  Department  of  Industrial 
Engineering  and  Operations  Research,  University  of  California, 
Berkeley. 
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EVERYMAN'S  GUIDE  TO  TIMES 


by 

Randall  R.  Willie 


TIMES  is  a system  of  computer  routines  for  time  series  analysis  and 
forecasting.  It  is  a tool  for  investigating  both  the  structure  of  a single 
time  series  and  the  relationship  between  a pair  of  series. 

The  system  is  based  on  a library  of  FORTRAN  subroutines  acquired  from 
the  Ohio  State  University  and  designed  for  analysis  of  time  series  data 
utilizing  techniques  of  Box  and  Jenkins  [1].  This  library  has  been  extended 
at  U.C.  Berkeley  to  include  a spectral  package  and  a user-oriented  control 
program.  In  addition,  the  original  subroutine  for  plotting  on  a line  printer 
has  been  replaced  by  more  sophisticated  routines  that  enhance  and  simplify 
graphical  output. 

The  control  program  accepts  time  series  data  and  a number  of  user  in- 
struction statements,  each  confined  to  a single  80-column  card.  The  purpose 
of  this  writeup  is  to  describe  the  control  program  input  options.  The  user 
should  have  some  familiarity  with  Box-Jenkins'  methods  and  terminology,  as 
presented  in  Reference  [1],  Examples  illustrating  program  output  are  pro- 
vided in  Appendices  II  and  III. 
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I.  GENERAL  INSTRUCTION  STATEMENT  FORMAT 

The  instruction  statements  all  have  a common  format.  Card  columns 
1 through  5 contain  a keyword,  which  is  a specific  string  of  five  or  less 
characters.  Columns  7 through  80  contain  an  entry  list,  which  consists  of 
numeric  character  strings,  or  entries,  separated  by  commas.  For  instance, 
this  card  is  a valid  instruction  statement: 

Col. 

} J 

IDEN  101,0,51 


The  keyword  is  'IDEN'  and  the  entries  are  '101',  'O',  and  '51'. 

Each  entry  represents  a numeric  value,  so  these  strings  may  consist 
only  of  the  characters  0 through  9,  plus  or  minus  signs,  a decimal  point, 
and  the  letter  E.  In  general,  any  numeric  character  string  acceptable  in 
the  FORTRAN  language  will  be  properly  interpreted  by  the  control  program. 

For  instance,  entries  that  may  be  used  are  integer  representations,  such 
as  '38.2'  or  '-.01',  or  scientific  notation,  such  as  '-1.0E-5'  '.IE-01', 
or  simply  'E25'.  Improper  strings  like  '1.0.0',  'FRED',  '+1.0+E01',  '.E01', 
'1.0E.1',  and  '+1E'  will  result  in  a SYNTAX  error.  Up  to  eight  blanks  may 
appear  immediately  preceding  or  following  an  entry  string  in  the  entry  list, 
but  blanks  may  not  appear  within  the  string;  for  example,  the  entries  '10  45', 
’ + 1.0E5',  or  '1.0  E5'  will  result  in  a SYNTAX  error. 

The  value  represented  by  each  entry  must  lie  between  an  upper  and  lower 
bound.  These  bounds  vary  depending  on  the  particular  keyword  and  entry;  a 
complete  list  is  provided  in  Appendix  I.  If  the  value  lies  outside  of  the 
bounds,  a BOUND  error  will  be  detected  by  the  control  program. 

Each  entry  value  is  assigned  to  a particular  parameter  in  the  control 
program,  and  the  control  program  determines  the  disposition  of  an  entry 


solely  by  the  number  of  commas  preceding  it  in  the  entry  list.  For  con- 
venience, many  of  the  entries  are  optional',  that  is,  the  user  may  choose 
to  leave  such  entries  out  of  the  list,  and  the  control  program  will  auto- 
matically assign  a default  value  to  the  parameters  associated  with  these 
entries.  However,  if  a particular  entry  is  left  out,  unless  all  of  the 
succeeding  entries  are  also  left  out,  the  comma  following  that  entry  must 
be  retained.  On  the  first  card  below,  the  first  entry  has  been  left  out; 
on  the  second  card,  the  second  entry: 

IDEN  ,0,51 
IDEN  101,, 51 

If  there  are  n possible  entries,  and  the  user  wishes  to  leave  out 
n - k (optional)  entries,  say,  the  portion  of  the  entry  list  after 
entry  may  be  left  out  entirely.  For  the  above  example,  if  the  last 
entries  are  ignored,  the  following  card  suffices: 

IDEN  101 

In  fact.  If  all  entries  are  optional,  an  instruction  card  need  only 
the  keyword: 

IDEN 

Unfortunately,  not  all  entries  for  every  keyword  are  optional.  Required 
entries  must  be  provided  by  the  user,  or  a REQUIRED  ENTRY  error  will  result. 

Output  from  the  control  program  includes  a listing  of  the  user  instruc- 
tion statements.  If  SYNTAX,  BOUND,  or  REQUIRED  ENTRY  errors  are  encountered, 
appropriate  messages  will  appear  below  the  instruction  statements  containing 
these  errors.  Execution  of  the  time  series  routines  will  not  begin  if  any 


the  last 

. i th 
the  k 

three 


contain 


of  these  errors  occur. 
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H.  MODELLING  AND  FORECASTING  WITH  A SINGLE  TIME  SERIES 

There  are  a number  of  approaches  to  extrapolation  of  a time  series, 
but  the  TIMES  program  is  primarily  suited  to  one  very  general  approach 
suggested  by  Box  and  Jenkins.  Based  on  the  first  and  second  order  sample 
moments  of  a given  series,  a particular  model  is  chosen  from  the  class  of 
autoregressive  integrated  moving  average  (ARIMA)  processes;  the  model  pa- 
rameters are  then  estimated  by  nonlinear  least  squares.  The  statistical 
adequacy  of  this  fitted  model  in  describing  the  data  may  be  judged  by  tech- 
niques discussed  in  Chapter  8 of  Reference  [1].  At  this  point  in  the  analy- 
sis, it  may  be  necessary  to  modify  the  form  of  the  model  and  reestimate  the 
parameters.  When  the  analyst  feels  that  a satisfactory  "fit"  to  the  data 
has  been  obtained,  the  model  can  be  used  for  forecasting.  The  principal 
stages  in  this  approach  are  model  identification } parameter  estimation , and 
forecasting. 

The  TIMES  control  program  permits  some  flexibility  in  performing  these 
stages.  The  control  program  processes  instruction  statements  and  time  series 
data  in  two  phases.  In  the  first  phase,  instruction  statements  are  decoded 
and  entry  values  and  time  series  data  are  stored.  The  second  phase  begins 
when  either  the  control  program  input  stream  ends  or  a RERUN  instruction  is 
encountered.  If  instruction  statement  errors  have  been  detected  in  the  first 
phase,  execution  terminates;  otherwise,  some  additional  error  checks  are 
made  to  insure  that  the  input  is  logically  consistent.  Finally,  the  time 
series  routines  are  executed.  If  the  input  stream  contains  an  IDEN  instruc- 
tion, the  identification  stage  is  performed  first.  The  estimation  stage  is 
performed  next  if  an  EST  instruction  has  been  provided.  The  FCST  instruction 
controls  the  forecasting  stage. 

One  execution  of  both  phases  of  the  control  program  will  be  called  a 
run.  Unless  the  RERUN  instruction  is  used,  a program  run  will  thus  be  the 
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same  as  a single  computer  job,  and  any  or  all  of  the  stages  of  Identification, 
estimation,  or  forecasting  may  be  performed.  The  order  of  Instruction  state- 
ments in  the  input  stream  for  a given  run  is  completely  arbitrary.  Instruc- 
tions having  certain  keywords  may  not  be  used  more  than  once  for  a run;  for 
example,  multiple  usage  of  the  keywords  IDEN,  EST,  and  FCST  is  not  permitted. 
Appendix  I specifies  for  each  keyword  whether  multiple  usage  is  allowed. 

Data  Input 

The  DATA  instruction  prepares  the  control  program  to  accept  data  values 
for  the  time  series.  The  format  of  this  instruction  is 

DATA  npts 

The  entry  npts  gives  the  number  of  data  points  in  the  series.  The  DATA  in- 
struction must  be  immediately  followed  by  a card  describing  the  data  format. 

A standard  FORTRAN  format  item  enclosed  in  parentheses  should  appear  on  this 
card,  beginning  in  column  1.  An  example  is 

(8F10.0) 

The  format  item  will  usually  be  of  the  form  aFb.O,  where  b denotes  the  width 
in  character  positions  of  a data  field,  and  a is  the  number  of  fields  on  a 
card. 

The  data  cards  follow  the  format  card.  Values  are  placed  successively 
in  data  fields  from  left  to  right  across  a card;  as  many  cards  may  be  used 
as  necessary  to  accommodate  all  time  series  points.  However,  all  cards  but 
the  last  must  have  a data  value  in  each  field.  This  value  should  contain  a 


decimal  point  and  may  appear  anywhere  within  the  field,  unless  E-notation  is 
used  (such  as  1.0E-02).  In  this  latter  case,  the  data  item  must  be  right- 
justified  in  its  field. 


r 


mm 


calculation  of  the  chi-square  statistic  to  test  for  white  noise. 

A single  program  run  for  the  model  identification  phase  may  be  made 
with  only  a DATA  instruction  and  data  package  and  an  IDEN  instruction.  The 
user  may  also  provide  an  SD1F  or  TRAN  card  if  seasonal  differencing  or  trans- 
formation of  the  input  series  ir  desired  prior  to  identification.  These  in- 
structions are  described  below. 


Differencing  and  Transformation 


Very  often  the  input  series  will  not  appear  stationary  or  will  contain 


strong  periodicities.  The  user  may  choose  to  apply  a transformation  or 
several  difference  operations  to  the  input  series  to  yield  a stationary 
series  that  can  be  suitably  modelled  as  an  autoregressive  moving  average 
(ARMA)  process.  More  specifically,  the  stationary  series  x^_  is  related 
to  the  input  series  in  the  following  manner: 


"s  / °i)ki 
x = n (i  - b f 

Li=1 


(1  - B)  ”T(Xt) 


Here  B is  the  familiar  backshift  operator  of  Reference  [1];  i.e.,  Bz  = ^ 

The  first  group  of  factors  in  the  above  expression  represent  seasonal  differ- 
ences, the  iC^  factor  denoting  k^  seasonal  differences  of  order  o^  , where 

o > 2 to  distinguish  seasonal  differences  from  regular  differences.  The 
k 

(1  - B)  ° factor  indicates  kQ  regular  differences.  Finally,  T(Xf)  is  a 
transformation  that  may  either  be  of  the  form  T(Xt)  = log(Xt  + t)  or 
T(X^)  * (Xt  + t)01  , 0 <_  a <_  1 . The  constant  t simply  insures  that 

X^  + r _>  0 , for  all  t . 

The  instruction  keywords  for  regular  differencing,  seasonal  differenc- 
ing, and  transformation  are,  respectively,  RDIF,  SDIF,  and  TRAN.  The  RDIF 


instruction  has  the  format: 


RDIF  [nrdiff] 


Nrdiff  is  the  number  of  regular  differences.  The  SDIF  instruction  is 


SD1F  ord,  [nsdiff] 


Ord  is  the  order  of  the  seasonal  difference  factor;  that  is,  the  integer 


Nsdiff  is  the  number  of  sea- 


quantity  o in  the  expression 


sonal  differences.  It  should  be  clear  that  no  more  than  one  RDIF  instruction 


will  be  necessary;  however,  as  many  as  four  SDIF  instructions  are  permitted 


The  TRAN  instruction  format  is 


TRAN  [alpha]  , [tau] 


Tau  is  any  quantity  t such  that  X + t _>  0 for  all  t . Alpha  Is  a value 


lying  in  [0,1]  . If  some  value  a > 0 is  chosen,  then  the  input  series 


will  be  transformed  according  to  T(X  ) = (X  + t) 


represent  the  initial  part  of  model  specification;  specification  is  completed 


by  choosing  a suitable  ARMA  model  for  the  x series.  The  RDIF  instruction 


is  not  effective  in  the  identification  stage,  but  SDIF  and  TRAN  instructions 


will  allow  the  user  to  identify  the  series  that  results  from  operating  on  the 


Model  Specification 


The  purpose  of  the  identification  stage  is  to  give  the  user  sufficient 


information  about  the  input  series  to  select  an  appropriate  ARIMA  model  which 


describes  the  process.  A set  of  instruction  cards  which  specify  the  form  of 


this  model  is  required  for  estimation  and  forecasting.  The  discussion  here 
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will  center  on  the  procedure  for  specifying  an  ARMA  model  for  the  stationary 
series  obtained  from  the  input  series  through  some  differencing  and 

transformation  scheme.  If  differencing  and  transformation  are  unnecessary, 
then  xt  will  denote  the  input  series  itself. 

The  TIMES  program  can  accommodate  a general  ARMA  model  of  the  form: 


II  <j>  (B)  (x  - u) 
i=l  1 C 


0 

6 + n 0 (B)a 

j=l  2 


Here  u is  the  mean  parameter  and  0^  is  the  trend  parameter.  Only  one 
of  these  two  parameters  will  appear  in  a particular  model;  if  xfc  is  the 
result  of  differencing  the  input  series,  then  no  mean  is  allowed  and  the 
trend  parameter  is  optional;  otherwise,  no  trend  is  allowed,  and  a mean 
parameter  is  required.  The  terms  ^(B)  and  0 (B)  represent,  respectively, 
the  1 autoregressive  factor  and  the  j moving  average  factor.  They  are 

polynomials  in  the  backshift  operator  B : 


/ ai  bi  Ci\ 

4>  (B)  *=  (1  - <p  B 1 - B 1 - ...  - $ B 1 ) 
i \ a.  bi  ci  / 

0 . (B)  =/l-0  B j - 0 Bj-...-0  BJ^ 

j V qj  rj  sj  / 


Finally,  a is  a pure  (white)  noise  series. 

A mean  or  trend  parameter  is  indicated  by  an  instruction  card  bearing 
the  keyword  MEAN  or  TRND.  Bolv>  cards  have  the  same  format: 


MEAN 

TRND 


[initv] 


For  the  estimation  stage,  initv  will  be  the  initial  parameter  value,  giving 


r 
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one  parameter  space  coordinate  location  where  the  nonlinear  least  squares 
minimization  search  is  to  begin.  If  the  forecasting  stage  is  performed  in 
the  same  program  run  as  estimation,  the  final  parameter  value  resulting  from 
the  minimization  process  will  be  passed  automatically  to  the  forecasting, 
subroutine.  However,  if  forecasting  is  done  in  a separate  run  from  esti- 
mation, initv  must  be  the  final  parameter  value  determined  from  the  esti- 
mation stage. 

The  keywords  AR  and  MA  are  used  for  autoregressive  and  moving  average 
parameters,  respectively.  A parameter  of  either  type  can  be  uniquely  iden- 
tified by  giving  the  index  of  the  factor  to  which  it  belongs  and  the  power 
of  B which  it  multiplies  in  that  factor.  The  indexing  scheme  for  the  fac- 
tors is  quite  arbitrary,  the  only  restrictions  being  that  each  factor  must 
have  an  integer  index  between  1 and  5,  and  no  two  factors  of  the  same  type 
(autoregressive  or  moving  average)  may  have  the  same  index.  Thus  up  to  5 
autoregressive  factors  and  5 moving  average  factors  may  appear  in  the  same 
model.  The  integer  power  of  B multiplied  is  the  entry  ovd  in  the  follow- 
ing format,  and  ifae  is  the  factor  index: 


i 


i 


it 


F 


SDIF  7 
TRAN  0 
AR  .1,1,1 
TRND  .2 
MA  .5,1,1 

MA  .8,7,2 


Estimation 

The  EST  instruction  informs  the  control  program  that  routines  for  pa- 
rameter estimation  are  to  be  executed.  The  format  is: 

EST  [nracor],  [rplot],  [nrchi],  [niter],  iepsl],  [eps2] 

After  estimates  of  model  parameters  have  been  found,  the  residual  series  is 
calculated  and  automatically  analyzed  by  the  same  routines  utilized  in  the 
identification  stage.  Thus  the  entries  nracor,  rplot , and  nrchi  have  the 
same  significance  with  respect  to  the  residual  series  as  the  entries  nacor, 
dplot , and  nchi  described  above.  Niter  gives  the  maximum  number  of  itera- 
tions permitted  in  the  nonlinear  (Marquardt)  minimization  routine.  Epsl 
is  the  relative  error  in  the  sum  of  squares  value  that  is  tolerable  when 
terminating  the  minimization  routine.  Similarly,  eps2  is  the  relative  error 
tolerable  in  the  estimated  value  of  each  model  parameter  when  terminating 
the  minimization  routine.  Output  from  the  estimation  phase  consists  of  a 
summary  of  the  status  of  each  iteration  of  the  minimization  routine,  fol- 
lowed by  the  estimated  parameter  values.  A cumulative  periodogram  and 
histogram  of  the  residuals  a^e  also  provided,  as  well  as  information  neces- 
sary for  the  identification  of  the  residuals. 

Forecasting 

The  forecasting  stage  is  controlled  by  a FCST  instruction: 

FCST  [conf  ] , nfcsts,  [to^,  [to2l,  [to^l,  [to^],  [to,.] 


Conf  may  be  one  of  the  integers  1,  2,  3,  4,  or  5,  which  indicate  that,  re- 
spectively, 50%,  75%,  90%,  95%,  or  99%  confidence  limits  are  to  be  provided 
with  the  forecast.  Nfosts  is  the  number  of  forecasts  to  be  made  forward 
from  each  of  the  time  origins  to^,  •••»  to^  • Though  the  entry  for  each 
time  origin  is  optional,  forecasts  are  done  only  for  consecutive  entries. 

For  instance,  the  following  card  will  result  in  forecasts  only  from  time 
origins  120  and  140: 

FCST  3,  15,  120,  140, ,160 

The  UPDA  instruction  permits  updated  forecasting  from  tne  kot  origin 
entered  on  the  FCST  card: 

UPDA  nupts,  [iuplot] 

This  instruction  may  only  be  used  once  in  a program  run  and  must  be  accom- 
panied by  a FCST  card.  Nupts  is  the  total  number  of  data  points  to  be  used 
in  updating;  it  should  be  noted  that  the  sum  of  the  last  time  origin  on  the 
FCST  card  and  nupts  may  not  exceed  the  number  of  data  points,  npts , in  the 
data  package.  For  the  i^  update  iteration,  points  x^,  ...,  xfcQ  +i  are 

assumed  to  be  known,  where  to^  is  the  last  time  origin.  Forecasts  are 
calculated  for  time  indices  to^  + i + l,to^  + i + 2,  ...,  to^  + nfcsts  . 
Forecast  values  and  confidence  limits  are  prir.ced  out  for  each  update  itera- 
tion; but  plotting  of  the  forecasts  is  optional,  since  it  will  usually  be 
expensive  and  unnecessary  to  produce  a new  graph  for  each  additional  point 
used  in  updating.  Iuplot  is  an  integer  between  1 and  nupts  that  controls 
the  number  of  forecast  plots  that  are  to  be  made.  If  iuplot  is  4,  for  ex- 
ample, a plot  is  made  after  each  group  of  four  successive  new  data  points 
has  been  incorporated  into  an  updated  forecast.  A plot  is  always  made  after 
all  nupts  have  been  utilized. 
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Multiple  Runs 

The  RERUN  instruction  card  has  no  entry  list  and  consists  only  of  the 
keywsrd  'RERUN'  in  columns  1 through  5.  The  main  purpose  of  this  card  is 
to  allow  the  estimation  and  forecasting  stages  to  be  performed  for  a number 
of  different  time  series  models  within  the  same  computer  job.  The  instruc- 
tion reinitializes  all  control  program  switches  and  storage  locations  except 
those  associated  with  the  DATA  and  NAME  instructions.  Thus,  following  a 
RERUN  instruction,  a new  set  of  model  specification  cards  may  be  provided 
as  well  as  EST  and  FCST  instructions;  and  the  time  series  routines  will  again 
be  executed,  utilizing  the  original  data  package  and  project  name.  If  a NAME 
or  DATA  instruction  is  used  following  the  RERUN  instruction,  the  new  name  or 
data  package  will  replace  the  original  information. 


* 

t 

3 


■ 
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III.  SINGLE  INPUT  TRANSFER  FUNCTION  MODELLING 

The  objective  of  transfer  function  modelling  is  to  investigate  the  re- 
lationship between  two  time  series.  It  is  assumed  that  the  dependent  series, 
denoted  here  by  y , can  be  represented  by  a linear  operation  on  the  inde- 
pendent series,  xt  , plus  a noise  component  nt  , where  y , , and  r_t 

are  all  stationary.  In  general,  xt  and  y^  will  be  generated  by  differ- 
encing and  possibly  transforming  the  actual  data  series  and  . Hope- 

fully, additional  information  available  from  the  Xfc  series  can  be  used  to 
obtain  better  forecasts  of  Y than  would  be  the  case  were  Y treated 
alone.  The  justification  and  motivation  for  the  techniques  considered  here 
are  developed  in  some  detail  in  Chapter  11  of  Reference  [1]. 

Transfer  function  modelling  involves  the  familiar  stages  of  identifica- 
tion, estimation,  and  forecasting.  These  transfer  function  stages  are  handled 
by  the  control  program  in  a manner  that  parallels  the  single-series  case. 

The  procedures  discussed  below  differ  from  those  of  Part  II  mainly  in  the 
complexity  of  model  specification. 

Data  Input  ■ 

The  DATA  and  NAME  instruction  of  Part  II  are  to  be  used  for  entering 
information  concerning  the  X series.  YDATA  and  YNAME  are  the  parallel 
instructions  associated  with  the  Y series;  these  instructions  have  the 
same  format  and  are  used  in  the  same  manner  as  DATA  and  NAME.  Should  the 
number  of  data  points  in  X differ  from  the  number  in  Yfc  , the  minimum 
of  the  two  series  lengths  will  be  used  for  analysis,  with  the  remaining 
points  in  the  longer  series  ignored. 

The  data  points  for  both  series  are  required  in  each  stage  of  analysis. 
The  NAME  and  YNAME  instructions  are  optional. 


16 


Differencing  and  Transformation 


The  scheme  of  regular  and  seasonal  differencing  and  data  transformation 
necessary  to  generate  the  stationary  series  xfc  and  y^  from  the  original 
series  Xt  and  Yfc  should  be  determined  prior  to  transfer  function  iden- 
tification. Differencing  and  transformation  for  Xt  can  be  determined  when 
this  series  is  analyzed  by  the  procedures  of  Part  II  to  obtain  the  prewhit- 
ening model,  which  is  discussed  below.  To  select  the  proper  differencing 
and  transformation  for  Y , it  may  be  desirable  to  use  Yfc  as  input  to  the 
single  series  identification  stage  of  Part  II. 

The  instructions  RDIF,  SDIF,  and  TRAN  introduced  earlier  operate  on  the 
Xt  series  in  transfer  function  modelling.  The  corresponding  instructions 
for  the  Y series  have  the  keywords  YRDIF,  YSDIF,  and  YTRAN,  and  the  for- 
mats of  these  cards  are  the  same,  respectively,  as  those  of  RDIF,  SDIF,  and 
TRAN.  Each  of  the  instructions  RDIF,  YRDIF,  TRAN,  and  YTRAN  may  be  used  only 
once.  However,  up  to  four  SDIF  instructions  and  four  YSDIF  instructions  are 
permitted  in  a run. 


Prewhitening  Model  Specification 


Transfer  function  modelling  involves  three  distinct  model  specifications: 
the  prewhitening  model,  the  noise  model,  and  the  transfer  function.  The  re- 
sults of  the  transfer  function  identification  stage  will  suggest  the  forms 
of  the  noise  model  and  the  transfer  function;  parameter  values  for  these 
two  models  will  then  be  determined  simultaneously  in  the  estimation  stage. 

The  prewhitening  model,  on  the  other  hand,  must  be  selected  prior  to  the 
identification  stage  by  considering  X alone  in  the  analysis  described  in 
Part  II.  The  residuals  from  applying  the  prewhitening  transformation  to  xt 
should  approximate  a pure  noise  series,  and  this  is  in  fact  the  objective  of 
single  series  analysis.  The  model  specification  discussed  in  Part  II, 


utilizing  the  keywords  AR,  MA,  MEAN,  and  TRND,  applies  without  modification 
to  the  prewhitening  model.  The  initv  entry  for  each  parameter  card  should 
be  the  final  estimated  value  obtained  for  that  parameter. 

The  same  prewhitening  model  instructions  should  be  provided  in  each 
stage  of  transfer  function  analysis.  If  no  prewhitening  model  instructions 
are  used,  the  prewhitening  transformation  will  be  the  identity. 

Identification 

The  stationary  series  x^  and  yt  will  result  from  differencing  and 
transforming  the  input  series  Xt  and  Yt  ; if  differencing  is  unnecessary, 
xt  an!  yt  will  be  the  original  or  transformed  series  translated  so  that 
the  means  are  zero.  (The  routines  perform  this  translation  automatically.) 

As  discussed  in  Reference  [1],  the  transfer  function  model  in  its  basic  form 
can  be  written  as 

yt  = Y + v(B)xt  + nt  , 

where  nf  is  a stationary  noise  series  with  zero  mean,  y is  a possible 
trend  parameter,  and  v(B)  is  an  operator  of  the  form 

b(B)  = |vq  + v^B  + v 2 B^  + ... 

The  purpose  of  the  identification  stage  is  to  provide  rough  estimates  of  the 
v weights  and  the  noise  series.  If  the  estimated  noise  series  has  a nonzero 
mean,  then  this  mean  can  be  taken  as  an  initial  estimate  of  y . 

When  the  prewhitening  transformation  is  applied  to  xt  , an  approximate 
pure  noise  series  will  result.  The  key  to  estimating  the  v weights 

is  to  apply  the  same  transformation  to  yfc  to  obtain  a residual  series  Bfc  ; 
it  is  shown  in  Reference  [1]  that  if  r Q(k)  is  the  sample  crosscorrelation 
function  between  a and  B„  , and  s and  s are  the  sample  variances. 
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raB(k)sB 


k * 0,1,2,  . . . 


An  estimate  of  the  noise  series  nt  may  then  be  calculated  according  to 
the  relation 

nt  + y = yt  - (vn  + vlB  + v2b2  + ...)xt  . 

The  TIDEN  instruction  card  initiates  the  above  manipulations  and  has 
the  format: 

TIDEN  [nccor],  [dplot],  [nvwts],  [nnacor],  [ndplot],  [nnchi] 


Nccor  is  the  number  of  lags  for  which  the  sample  crosscorrelation  function 

r (•)  is  to  be  calculated.  If  dplot  is  1,  a combined  plot  of  the  original 
etc 

series  Xt  and  Y is  provided.  The  entry  nvwts  specifies  the  number  of 
v-weight  estimates  to  be  utilized  in  calculating  the  noise  series.  The  es- 
timated noise  series  is  analyzed  by  the  single  series  identification  routine 
of  Part  II.  Thus  nnacor , ndplot , and  nnchi  correspond,  respectively,  to  the 
entries  nacor , dplot,  and  nchi  on  the  IDEN  instruction  card. 


Transfer  Function  Specification 

The  v weights  are  used  only  indirectly  in  formulating  the  transfer 
function  model;  in  the  simplest  case,  these  weights  are  taken  to  be  functions 
of  a small  number  of  parameters  associated  with  an  input  lag  factor,  oo(B)  , 
and  an  output  lag  factor,  6(B)  . These  factors  are  of  the  form 

w(E)  = ( u)_  - uj,  B - ...  - „ BP 
\ 0 1 p 

6(B)  = (l  - 6lB  - ...  - 6qBqj  , 
and  v(B)  is  assumed  to  have  the  representation 
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v(B)  = 5-1(B)w(B)Bb 

Here  b will  be  greater  than  zero  If  = • . . = ^ = 0 . In  general, 

oj(B)  and  6(B)  are  to  be  chosen  with  as  few  parameters  as  will  allow  the 
coefficients  of  B^  in  6 ‘*‘(B)<u(B)Bb  to  approximate  the  estimated  sequence 
determined  in  the  identification  stage.  The  process  of  selecting  input 
and  output  lag  factors  may  require  considerable  judgment  and  care  on  the  part 
of  the  analyst.  This  process  will  also  yield  initial  estimates  of  parameters 

, i = 1,  ...,  p and  6 , j = 1 c . 

The  estimation  and  forecasting  routines,  in  fact,  can  accommodate  mul- 
tiple lag  factors,  since  seasonal  ’relationships  between  y^  and  may 

require  this  additional  flexibility.  Thus  the  general  model  is  as  follows: 

n . n 

6 a) 

y = Y + n 6 (B)  n OJ  (B)x  , + n , 

c j=l  3 i=l  1 t_D 

where 


^(B) 

m^(B) 

VB> 


/ *i  bi 

(1-U)  B - ui  B - 

l *i  hi 

/ ql  rl 

I a).  - 0)  B - u B 

l °1  ql  rl 

(l  - 6 BUJ  - 6 B^  - 

V uj  vj 


w B , i ^ 1 
Ci  / 

- w B 
s. 


The  number  of  factors  of  either  type  is  limited  to  five.  Each  transfer 
function  parameter  can  be  uniquely  identified  by  the  index  of  the  lag  factor 
to  which  it  belongs  and  the  power  of  B which  it  multiplies.  It  should  be 
noted  that  the  only  parameter  that  does  not  multiply  a positive  power  of  B 


i 
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is  u)  ; that  is,  this  parameter  must  be  in  the  first  input  lag  factor. 

U1 

With  this  exception,  the  indexing  of  the  lag  factors  is  unimportant. 

The  instruction  cards  for  the  input  lag  and  output  lag  parameters  have 

the  respective  keywords  ILAG  and  OLAG: 

I LAG 
OLAG 


[initv],  ord,  ifac 


Initv  is  the  initial  value  of  the  parameter  for  estimation.  As  before,  if 
forecasting  is  done  in  a separate  run  from  estimation,  initv  must  be  the 
value  of  the  parameter  from  estimation.  Ord  is  the  power  of  B that  the 
parameter  multiplies,  and  ifac  is  a number  between  1 and  5 giving  the  index 
of  the  factor  to  which  the  parameter  belongs. 

Transfer  function  model  specification  is  completed  by  optionally  speci- 
fying the  trend  parameter  y and  the  number  of  lags  b by  which  perturba- 
tions in  y lag  those  in  xt  . These  two  instructions  are  as  follows: 

TTRND  [initv] 

B b 

Up  to  ten  B instructions  are  permitted  in  an  estimation  run.  Parameter 
estimation  will  be  repeated  for  each  b value. 

Noise  Model  Specification 

The  noise  series  nt  is  assumed  to  be  stationary  and  have  zero  mean, 
but  in  most  cases  it  will  not  be  a pure  noise  series.  An  analysis  of  the 
sample  autocorrelation  and  partial  autocorrelation  functions  of  the  esti- 
mated series  nt  calculated  in  the  identification  stage  will  suggest  an 
appropriate  model  for  n^  in  terms  of  autoregressive  and  moving  average 
factors.  The  general  form  of  the  model  is 


n <MB)n  = n 0 (B)a  , 
i-i  c j«i  3 z 
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where  at  is  a pure  noise  series  and 


/ 3i  bi  Ci\ 

.(B)  - [1  - * B A - ♦ B 1 - ...  - <t>  B 1 

1 \ ai  bi  Ci  / 

. (B)  = (l  - 0 B^-e  B J - ...  - 6 B^l 

J l qj  rJ  SJ  ) 


Thus  each  parameter  in  the  noise  model  is  either  an  autoregressive  or  a 
moving  average  parameter  and  is  uniquely  identified  by  the  index  of  the 
factor  to  which  it  belongs  and  the  power  of  B which  it  multiplies.  The 
keywords  are  NAR  and  NMA  for  autoregressive  and  moving  average  parameters, 
respectively,  and  the  instruction  card  formats  are 


NAR 

NMA 


[initv],  ord,  ifac 


InitVj  ordj  and  ifac  are  defined  in  the  usual  manner,  and  up  to  five  factors 
of  either  type  are  permitted. 

A Specification  Example 

An  example  is  appropriate  to  illustrate  the  use  of  model  specification 
instructions.  The  structure  of  the  two  series  considered  in  Section  11.5.3 
of  Reference  [1]  may  be  expressed  as 

y = .035  + (1  - .72B)-1(4.82)xt_3  + (1  - .54B)at 


xt  = (1  - .32B)at  , 


where  and  xt  are  first  differences  of  the  original  series  and  a^ 

and  c*t  are  pure  noise  series.  The  complete  set  of  instruction  cards  for 
all  three  models,  plus  differencing  instructions,  is  as  follows: 
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RDIF 

YRDIF 

MA 

.32,1,1 

OLAG 

•72,1,1 

I LAG 

A. 82, 0,1 

B 

3 

NMA 

•54,1,1 

Estimation 

The  estimation  stage  is  very  similar  to  the  single  series  estimation 
of  Part  II.  The  prewhitening,  transfer  function,  and  noise  models  should 
be  specified,  with  initial  estimates  for  transfer  function  and  noise  para- 
meters. The  instruction  format  for  estimation  is 

TEST  [nrccor],  [nracor],  [rdplot],  [nrchi],  [niter],  [epsl],  [eps2] 

The  first  five  entries  are  associated  with  the  analysis  of  the  residuals, 
at  , from  the  estimated  noise  series  nt  , which  is  determined  after  final 
parameter  estimates  have  been  obtained.  Nracor  is  the  desired  number  of 
lags  for  which  the  sample  crosscorrelation  function  between  at  and  the 
prewhitened  xt  series,  , is  calculated.  As  explained  in  Section  11.3.3 
of  Reference  [1],  this  crosscorrelation  is  a measure  of  the  discrepancy  be- 
tween the  "correct"  v weights  for  the  transfer  function  model  and  the  v 
weights  that  are  calculated  from  the  parameter  estimates.  The  transfer 
function  estimation  routine  utilizes  the  single  series  identification  rou- 
tine to  provide  information  for  analysis  of  the  residuals;  so  again,  nracor, 
rdplot,  and  nrchi  correspond  to  the  entries  for  the  IDEN  instruction. 

Finally,  niter,  epsl,  and  eps2  are  control  values  for  terminating  the 
le  st-squares  minimization  process;  they  are  considered  in  Part  II  in 
conjunction  with  the  EST  instruction. 

Forecasting 


For  the  forecasting  stage,  all  three  models  should  be  specified,  and 
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if  this  stage  is  run  separately  from  estimation,  the  initv  entry  for  each 
transfer  function  and  noise  parameter  must  be  the  final  estimate  for  that 
parameter.  TFCST  is  the  forecasting  instruction: 

TFCST  [conf  ] , [nfcsts],  [to^,  [to,,],  [to3],  [to^],  [to5] 

The  entries  for  this  instruction  correspond  exaccly  to  the  entries  on  the 
FCST  instruction  card.  The  single-series  forecasting  routine  automatically 
provides  forecasts  of  the  series,  utilizing  the  prewhitening  model. 

The  forecasts  are  then  taken  as  input  to  the  transfer  function  model 

to  eventually  yield  forecasts  of  Y . 

The  single-series  option  of  updated  forecasting  is  not  available  in 
the  transfer  function  forecasting  stage. 

Multiple  Runs 

The  RERUN  instruction  may  also  be  used  in  transfer  function  modelling. 
The  instruction  reinitializes  the  control  program,  but  leaves  data  and  names 
for  both  series  unchanged. 


i 
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IV.  SPECTPAL  IDENTIFICATION 

Spectral  methods  of  time  series  analysis  differ  significantly  from  the 
Box-Jenkins  approach,  but  these  methods  are  often  quite  useful,  especially 
for  time  series  arising  from  natural  phenomena,  such  as  seismic,  atmospheric, 
or  biological  processes.  Detailed  discussions  of  spectral  techniques  may  be 
found  in  Brillinger  [2]  and  Watts  and  Jenkins  [3].  The  TIMES  program  has  op- 
tions for  both  the  spectral  analvs  s of  a single  time  series  and  the  cross- 
spectral  analysis  of  a pair  of  series. 

Let  c (•)  be  the  autocovariance  function  of  a time  series  . The 
xx  t 

power  spectrum  of  this  series  is  the  Fourier  transform  of  the  autocovariance 
function: 

f (w)  = ■£-  £ c (k)e  0 < w < .5  , i = S-l 

XX  2 TT  £■  w ' ’ — — ’ 


k=-° 


XX 


Over  any  finite  time  interval,  the  process  xt  may  be  approximated  by  a 
process  x^  consisting  of  a sum  of  sinusoidal  components  having  random 
amplitudes  and  phases.  These  components  are  uncorrelated  and  have  frequen- 
cies in  the  range  -.5  to  .5  cycles  per  time  unit.  (A  time  unit  is  the  dura- 
tion between  successive  points  of  the  xt  series.)  The  approximation  can 

« 2 

be  made  arbitrarily  close  in  the  sense  that  E(xt  - xt)  can  be  made  small. 

„ 2 

As  E(xt  - xt>  0 , these  components  may  be  thought  of  as  a decomposition 
of  the  xt  process,  and  the  power  spectrum  represents  the  squared  component 
magnitudes  as  a function  of  the  frequency  w . 

In  Part  III  the  dependent  series  y was  assumed  to  be  related  to  xt 
by  > 

yt  “ y + v(B)xt  + nt  , 


w 
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where  y is  a possible  trend  parameter,  v(B)  is  the  generati: g function 

of  the  v-weights,  and  n£  is  a stationary  noise  series,  uncorrelated  with 

x . Let  c (•)  be  the  cross  covariance  function  between  y^  and  x_  ; 
t yx  t t 

that  is,  c^(k)  = CovIy^^.x^J  ' T^e  aross  spe°t'nirn  is  the  Fourier  trans- 
form of  the  cross  covariance: 


f (w)  * ~ £ c (k)e  ^^7T^CW  - .5  < w < .5  . 

W 9 IT  L \TY  


yx  2ir  ^ yx' 


It  is  not  difficult  to  show  that 


OO  f 

-12ttwn  r -i2rrkw  yx^ 


, -lzirwx  r 
v(e  ) = 1 vke 

k=-°° 


f (w) 

XX 


Now  v(e  *^1TW)  may  be  written  as  the  product  of  a real  function  G(w)  , 
called  the  gain  function,  and  e*2n<Kw)  ^ where  <},(w)  ±s  a (real)  phase 


function: 


, -i2iTO\  n . v i27nf>(w) 
v(e  ) = G(w)e 


The  sinusoidal  components  associated  with  v(B)xt  can  essentially  be  ob- 
tained by  multiplying  the  components  of  xt  by  v(e  ^^irw)  . This  multi- 
plication alters  the  magnitude  of  the  component  with  frequency  w through 

the  factor  G(w)  ; the  phase  is  shifted  by  the  factor  e*2ir<f>(w)  . The 

*“  i2  TTW 

pulse  response  function  may  be  evaluated  from  v(e  ) through  an  inverse 
Fourier  transform: 


'f 

I . -12 ttw,  i2ffwk, 

I v(e  )e  dw 


Another  important  function  in  cross-spectral  analysis  is  tne  coherence , 


defined  by 
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2 ,fyx(w)l 

*yx(w)  " f„_(v)f__.(w) 


XX 


yy 


Roughly,  this  quantity  behaves  like  a squared  correlation  coefficient  for 
each  frequency  w , measuring  the  correlation  between  the  sinusoidal  com- 


ponent of  yt  and  that  of  xt  at  the  frequency  w . Furthermore,  it  can 
be  shown  that  the  spectrum  of  the  noise  series  may  be  evaluated  as 


f (w)  = f (w) 
nn  yy 


[1  - 4<w>] 


Thus  ICyX(w)  is  in  some  sense  a measure  of  the  portion  of  information  in 
the  y series  that  is  directly  attributable  to  the  x^  series.  Finally, 
the  noise  auto covariance  function  may  be  obtained  by  taking  the  Fourier  in- 
verse of  f (•)  : 
nn 

1/2 

c (k)  -2  I f (w)cos(2wk)dw  . 
nn  J nn 

0 


For  a single  series  x^  , the  TIMES  spectral  analysis  routine  simply 

estimates  the  power  spectrum  f^C*)  • a pair  of  series,  xt  and  y , 

the  routine  first  estimates  f (•)  , f (•)  , and  f (•)  . These  estimates 

xx  yy  yx 

are  used  in  turn  to  estimate  the  gain,  phase,  coherence,  and  noise  spectrum, 

in  a manner  suggested  by  the  above  definitions  of  these  functions  in  terms 

of  f (•)  , f (•)  , and  f (•)  . From  the  gain  and  phase  estimates  the 
xx  yy  yx 

impulse  response  weights  are  determined  by  Fourier  inversion,  and  the  noise 

spectrum  estimate  is  inverted  to  yield  the  noise  autocorrelation  function. 

Estimates  of  f (■)  , f (•)  , and  f (•)  are  based  on  the  discrete 
xx  * yy  yx 

Fourier  transform defined  for  an  arbitrary  (real  or  complex)  time  series 


, t = 0,  • • . , T — 1 , by 
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dT(w)  - £ z e i2lttw  f - .5  < w < .5 
Z t-0  c " " 


The  Fast  Fourier  Transform  algorithm  is  used  to  compute  the  sample  discrete 

T N 

transform,  so  a particular  realization  of  dz(*)  is  calculated  only  for  2 

N N 

equally  spaced  frequencies  w = j/2  - .5  , j = 1,2,  ...,  2 , where  N is 
chosen  such  that  2N  ^ T . For  the  series  xfc  and  yt  , both  of  length  T , 
periodograms  1^(*)  and  I (•)  are  defined  by 


Txx(w)  = 2^T  dx(w)dx(w)  ’ 


Zyy(W)  = 2wT  dy(w)dy(w)  ’ 0 < w < .5  , 


where  the  bar  denotes  conjugation.  The  cross  periodogram  is 


Zyx(w)  = 2*T  dy (w)dx(w)  * ".5  1 w < .5 


Now  I^Cw)  » I (w)  » aad  I^Cw)  are  ran<*om  variables  that  have  respec- 
tive expected  values  f (w)  , f (w)  , and  f (w)  , as  shown  in  Reference 
r xx  yy  yx 

[2];  however,  these  periodogram  ordinates  have  variances  that  do  not  decrease 

T T T 

as  T increases.  More  stable  estimators  f (w)  , f (w)  , and  f (w)  are 

xx  yy  yx 

therefore  obtained  by  smoothing  the  periodograms: 


‘ <2mVl)  ” W”  + t/2N)  • 

m 


J_m  V” + k/2N)  • 


] 


4(»>  ■ 7STTT7  . l V"  + t/2N) 

k=-m 


The  quantity  (2m  + l)/2  is  the  smoothing  bandwidth , and  it  must  be  chosen 
large  enough  by  the  analyst  so  that  spectral  estimates  are  stable,  but  not 
so  large  as  to  obscure  details  of  the  spectra.  It  is  good  practice  to  per- 
form the  analysis  several  times  using  different  bandwidths,  most  commonly  on 
the  range  .01  to  .1. 

Reference  [2]  discusses  the  asymptotic  distributions  of  the  estimators 
T T T 

f (•)  , f (*)  and  f (•)  and  then  derives  various  distributional  proper- 
xx  yy  yx 

ties  for  the  coherence,  phase,  gain,  and  residual  spectrum  estimators.  Ap- 
proximate confidence  limits  for  the  estimates  can  therefore  be  calculated, 
and  these  approximations  should  be  quite  valid  for  data  samples  consisting 
of  more  than  about  250  points.  Estimates  of  power  spectra  and  gain  function 
are  more  appropriately  given  in  terms  of  logarithms,  since  this  simplifies 
computation  of  confidence  limits.  The  TIMES  spectral  routine  outputs  l°8^g 
spectra  and  gain. 

The  spectral  analysis  subroutine  is  called  if  an  SIDEN  instruction  is 
encountered  in  the  control  program  input  stream.  The  instruction  format  is 


SIDEN  [opt],  [bw..],  [bw,],  [bw.],  [bw,],  [bw  ] 


The  entry  opt  determines  the  number  of  frequencies  for  which  spectra,  gain, 
phase,  and  coherence  will  be  estimated.  This  entry  does  not  determine  the 
frequency  range,  which  is  always  [ 0, . 5 ] ; rather,  it  determines  how  finely 
this  interval  is  to  be  divided  for  discrete  output.  If  opt  is  1,  calcula- 
tions are  done  for  101  frequencies.  Opt  may  be  2 only  if  the  data  series 
passed  to  the  spectral  routine  exceed  512  poii ts  in  length,  and  in  this  case 
calculations  are  done  for  251  frequencies.  The  entries  bw i,  i - 1,  ...,  5 , 


• • • t 
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are  smoothing  bandwidths.  A complete  spectral  analysis  will  be  done  for 
each  successive  bandwidth  entry. 

Use  of  the  SIDEN  instruction  will  result  in  estimation  of  one  power 
spectrum  if  the  input  stream  contains  only  a single  data  package.  If  data 
packages  for  both  the  Xt  and  Y series  are  present,  cross-spectral 
analysis  will  be  done.  The  RDIF,  SDIF,  TRAN,  YRDIF,  YSDIF,  and  YTRAN  in- 
structions will  apply  differencing  and  transformation  operations  to  the 
series  prior  to  spectral  analysis.  In  addition,  a prewhitening  model  may 
be  specified  in  the  same  manner  as  discussed  in  Part  III;  the  prewhitening 
technique  should,  in  fact,  always  be  used  when  cross-spectral  estimates  are 
desired.  Single-series  spectral  analysis  gives  the  sample  power  spectrum 
of  the  residuals  from  applying  the  prewhitening  model  to  the  Xt  series, 
so  in  single-series  modelling  the  power  spectrum  is  an  additional  tool  to 
aid  the  analyst  in  identif ication  of  the  residuals.  If  an  SIDEN  and  an  FST 
instruction  are  used  in  the  same  run,  spectral  analysis  will  follow  the  esti- 
mation stage  in  the  processing  sequence;  thus,  final  parameter  estimates, 
rather  than  the  initial  values  on  AR  and  MA  instructions,  will  be  used  for 
the  prewhitening  transformation. 


APPENDIX  I 

Instruction  Entry  Bounds  and  Defaults 

The  quantities  NOBMX,  NACMX,  and  NFMX  are,  respectively,  the  maximum 
number  of  observations,  autocorrelations,  and  forecasts  that  are  permitted 
in  a particular  implementation  of  the  TIMES  program.  These  quantities  depend 
on  the  size  of  storage  arrays  in  the  control  program;  sizes  may  be  altered 
only  by  changing  several  FORTRAN  instructions  in  the  program. 

The  abbreviation  "req."  indicates  that  the  associated  entry  is  required 
and  no  default  value  can  be  assigned.  The  abbreviation  "succ."  indicates 
that  forecasting  or  spectral  analysis  will  be  done  only  for  successive  time 
origin  entries  or  bandwidth  entries.  The  "Multiple  Usage"  column  contains 
a "yes"  if  two  or  more  instructions  having  the  associated  keyword  are  per- 
mitted in  a program  run. 

It  should  be  noted  that  a run  may  not  be  successful  even  if  all  entries 
satisfy  the  bound  constraints  below:  The  control  program  offers  only  a pre- 


liminary  check  on  the  validity  of  some 

entries;  final 

checks  are 

done  when 

the  time 

series  routines 

are  executed. 

Multiple 

Lower 

Upper 

Keyword 

Usage 

Entries 

Bound 

Default 

Bound 

DATA 

no 

npts 

1 

(req.) 

NOBMX 

NAME 

no 

(none) 

IDEN 

no 

nacor 

1 

npts/5 

NACMX 

dplot 

0 

0 

1 

nchi 

1 

nacor 

NACMX 

RDIF 

no 

nrdiff 

1 

1 

100 

SDIF 

yes 

ord 

2 

(req.) 

io30 

nsdif f 

1 

1 

100 

TRAN 

no 

alph.' 

0 

0 

1 

tau 

-io30 

0 

io30 

I 

I 

J 


!»*«Q 
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Keyword 

Multiple 

Usatse 

Entries 

Lower 

Bound 

Default 

Upper 

Bound 

MEAN 

no 

initv 

-1030 

.1 

io30 

TRND 

no 

initv 

-io30 

.1 

io30 

AR 

yes 

initv 

-1 

.1 

1 

ord 

1 

(req. ) 

io30 

ifac 

1 

(req. ) 

5 

MA 

yes 

initv 

-1 

.1 

1 

ord 

1 

(req. ) 

100 

ifac 

1 

(req.) 

5 

EST 

no 

nracor 

1 

npts/5 

NACMX 

rplot 

0 

0 

1 

nrchi 

1 

nracor 

NACMX 

niter 

0 

15 

999 

epsl 

io-8 

icf 3 

.1 

eps2 

icf8 

io-3 

.1 

FCST 

no 

conf 

1 

3 

5 

nfcsts 

1 

(req.) 

NFMX 

tot  (1  1 i 1 5) 

1 

(succ. ) 

npts 

UPDA 

no 

nupts 

1 

(req. ) 

NFMX 

iuplot 

1 

nupts 

nupts 

RERUN 

(none) 

YDATA 

no 

npts 

1 

(req.) 

NOBMX 

YNAME 

no 

(none) 

YRDIF 

no 

nrdif f 

1 

1 

100 

YSDIF 

yes 

ord 

2 

(req. ) 

io30 

nsdif f 

1 

1 

100 

YTRAN 

no 

alpha 

0 

0 

1 

tau 

O 

CO 

O 

1 

0 

io30 

-----  ■■  ■ 


Ml  - - I - 
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Keyword 

Multiple 

Usage 

Entries 

Lower 

Bound 

Default 

Upper 

Bound 

TIDEN 

no 

nccor 

1 

npts/5 

NACMX 

dplot 

0 

0 

1 

nvwts 

1 

nccor 

NACMX 

nnacor 

1 

npts/5 

NACMX 

ndplot 

0 

0 

1 

nnchi 

1 

nnacor 

NACMX 

I LAG 

yes 

initv 

-1030 

.1 

io30 

ord 

0 

(req. ) 

1030 

ifac 

1 

(req. ) 

5 

OLAG 

yes 

initv 

-1 

.1 

1 

ord 

1 

(req.) 

io30 

ifac 

1 

(req.) 

5 

TTRND 

no 

initv 

-1030 

.1 

io30 

B 

yes 

initv 

-1030 

(req.) 

io30 

NAR 

yes 

initv 

-1 

.1 

1 

ord 

1 

(req.) 

1030 

ifac 

1 

(req.) 

5 

NMA 

yes 

initv 

-1 

.1 

1 

ord 

1 

(req.) 

io30 

ifac 

1 

(req.) 

5 

TEST 

no 

nrccor 

1 

npts/5 

NACMX 

nracor 

1 

npts/5 

NACMX 

r dp lot 

0 

0 

1 

nrchi 

1 

nracor 

NACMX 

niter 

0 

15 

999 

epsl 

10"8 

IQ'3 

.1 

eps2 

10-8 

lO'3 

.1 

j 


33 


Keyword 

Multiple 

Usage 

Entries 

Lower 

Bound 

Default 

Upper 

Bound 

TFCST 

no 

conf 

1 

3 

5 

nfcsts 

1 

(req. ) 

NFMX 

to.^  (1  <_  i <_  5) 

1 

(succ. ) 

npts 

SIDEN 

opt 

1 

1 

2 

bw±  (1  1 i 1 5) 

0 

(succ. ) 

.25 

APPENDIX  II 

An  Example  of  Modelling  and  Forecasting  with 
a Single  Time  Series 

Data  for  the  following  example  is  given  in  Reference  [1],  p.  531.  The 
model  for  this  particular  series  is  considered  in  some  detail  in  Chapter  9 
j of  Reference  [1]. 


INTER.  NATION  A U A1R.LINE  RASSENGEfilS 


450 


an* 


150 


soo 


ITERAtiON  NC 


♦ 

♦ 

I 

J 

j 

1 

| 

! 

I 

* 

I 

j 

♦ 

* 

i 

i 

♦ j 

i 

i 

j 

- 

♦ 

I 

! 

♦ i 

! 

1 

1 

♦ 

♦ 

| 

* I 

j 

♦ 

♦ 

1 

} 

| 

* « 

♦ 

1 

♦ 

• 

j 

I 

| 

• 

• 

• 

j 

• mm 

• 

• • m 

! 

• • « 

j 

• • • 

• • 

1 

• 

• •• 

♦ 

♦ i 

• • 1 

. -* 

« • « 

• t 1 

* 

| 

• » • 

• • • 

• 

« • < 

• •*  • 

• • • 

• • 

j 

• • m 

• • 4 

• 

* 

• «i 

• • • 

• > « 

• • « 

• • « 

•* 

o.  • • 

• « • 

• mm 

(%• 

a • ♦ 

• • H 

• • 

♦ 

a.  • « 

• • • 

• • « 

• • 

• a < 

• • 1 

• • 1 

• • a. 

• • 

• • • 

• 

• • « 

cl  • • 

• • ■ 

• • 

* 

• • • 

• • a 

• • 4 

• • « 

* 

• • « 

a • 

• • i 

*-• 

• • • 

• • H 

a • • 

• • 

> • < 

• • « 

• a « 

» • . « 

< « < 

M « 

•<Q.i 

4 • 

•»  • * 

* « ♦ 

• • . 

• ♦a 

* * 

• * 3 

• • CL 

• • • 

• • •• 

• • • 

• • ■ 

• • 1 

e • • 

• • « 

# • 

• . . 

• • • 

•L  J 

• • H 

* - 

• • • 

• • 4 

• * a. 

• • « 

>* 

« # 

•I  •*  < 

i •«  • a 

1 ml  4 4> 

-• 

• • • 

• • « 

• • Q. 

• • • 

• • * 

• • * 

• • 0. 

• • 4 

• • • 

• • 4 

• • 

a • * 

• • 

• • m 

• • 

n • « 

♦ 

■ • • 

• • 

♦ a • 

* 

« • • 

• • 

• a •. 

• 

• • 4 

• • 

• * ai 

• • 

• • CL 

• • • 

• • 

• • Qj 

••  • 

• • 

• • ^ 

• « 

• • 

• • «i 

• 

« • 

• • 4 

• 

• • 

• • 4 

• • 

• • 4 

♦ 

♦ 

• • 

4 • • 

«o 

>o 

*© 

*c 

•c 

•c 

r) 

CM 

f 

1 • 

« ♦ 

1 

t 

i 

STIH 


600 


300 


600 


001 


•NO  ipdatino  forecasts  AFTER  NEW  OBSERVATION 


59 


APPENDIX  III 

An  Example  of  Transfer  Function  Modelling 
and  Forecasting 

Data  for  the  following  example  is  given  in  Reference  [1],  p.  537.  A 
transfer  function  model  for  this  data  is  presented  in  Section  11.5.3  of 
Reference  [1]. 
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